stability of ultracold atomic Bose condensates with Rashba spin-orbit coupling against 

quantum and thermal fluctuations 
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We study the stability of Bose condensates with Rashba-Dresselhaus spin-orbit coupling in three dimen- 
sions against quantum and thermal fluctuations. The ground state depletion of the plane-wave condensate 
due to quantum fluctuations is, as we show, finite, and therefore the condensate is stable. We also calculate 
the corresponding shift of the ground state energy. Although the system cannot condense in the absence 
of interparticle interactions, we show by estimating the number of excited particles that interactions sta- 
bilize the condensate even at non-zero temperature. Unlike in the usual Bose gas, the normal phase is 
not kinematically forbidden at any temperature; calculating the free energy of the normal phase at finite 
temperature, and comparing with the free energy of the condensed state, we infer that generally the system 
is condensed at zero temperature, and undergoes a transition to normal at non-zero temperature. 



Creation of artificial gauge fields in ultracold atoms 
presents opportunities to study physical phenomena pre- 
viously unattainable in atomic or condensed matter sys- 
tems |1]. Of current interest is the non-Abelian Rashba- 
Dresselhaus spin-orbit coupling L2i ^ which, following 
several proposals MW, has recently been realized exper- 
imentally [8, 9]. Depending on the details of the Rashba- 
Dresselhaus coupling and interparticle interactions, one of 
two possible ground state phases is expected: a plane-wave 
state, in which the system condenses into a single momen- 
tum state, and a striped phase, where the condensate is a 
linear superposition of two states with opposite momenta 
[10, 11, reviewed in 12]; renormalization of the interpar- 
ticle interaction beyond mean-field 113, .14.1 has significant 
effects on the phase diagram 11131 flSll . 



Condensates of ultracold bosons in three dimensions 
with Rashba spin-orbit coupling differ from usual Bose- 
Einstein condensates (EEC's) in several important ways. 
In the absence of interparticle interactions, the low-lying 
density of states is two-dimensional iflql . and thus conden- 
sation is destroyed by thermal fluctuations at any non-zero 
temperature. With interparticle interactions present, fluc- 
tuations around mean-field states lead at finite temperature 
to an instability of the plane-wave state in two dimensions 
iflTll . In this paper, we focus on three-dimensional ultracold 
bosons with isotropic Rashba-Dresselhaus coupling in the 
x-y plane, to investigate the effects of quantum and thermal 
fluctuations on a plane-wave Bose-Einstein condensation, 
and show that interactions in fact stabilize the condensate 
in 3D. This interaction-induced EEC is a unique feature of 
bosons with Rashba-Dresselhaus spin-orbit coupling, with 
no analogous system yet found. However, unlike in usual 
EEC's, a non-condensed state is not, as we show, kinemat- 
ically forbidden at any non-zero temperature. Condensa- 
tion, while favored at very low temperature, should dis- 
appear at high temperature. As in a ECS superconduc- 
tor, where both a normal and condensed state are allowed 
at low temperature, the system should undergo a similar 
phase transition at a critical temperature. 



We consider bosons with an isotropic Rashba- 
Dresselhaus spin-orbit coupling in three dimensions 
with an isotropic interaction, described by the Hamiltonian 
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As previously lll4l [iSll . m is the atomic mass, V is the 
volume of the system, and k is the spin-orbit coupling 
strength, taken to be positive. The isotropic s-wave cou- 
pling is g. The operators Op and 6p annihilate atoms with 
momentum p in the pseudospin states a and 6, respectively. 
The a^ and a.y are the usual Pauli matrices between the in- 
ternal states, and / is the two-by-two identity matrix. The 
dispersion relation of the single particle terms in the Hami- 
tonian has two branches e±(p) = {{p±_ ± kY +pl}/2m, 

where p^ = ^/p1 + pi, with circularly degenerate ground 

states along {p_L,Pz) = ('^, 0). In this paper, we assume 
that g is the (constant) mean-field coupling; extension be- 
yond mean-field coupling, as in il3Ml5h . is left as a future 
problem. Starting from the plane -wave ground state with 
momentum k, = (k, 0, 0), we construct the single-partice 
Green's functions, including quantum fluctuations via the 
Eogoliubov approximation, and derive the low-momentum 
spectra, estimate the condensate depletion, and calculate 
the ground state energy as a function of the spin-orbit cou- 
pling strength. As we show, the number of excited parti- 
cles increases and the ground state energy decreases with 
increasing spin-orbit coupling strength. 

Since the operator (a], — 5], )/v2 creates a particle in the 
single-particle ground state with momentum k, it is easier 
to work in the following (— , +) basis: 
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The state created by V'-.k is macroscopically occupied. 



We first derive the fluctuations of the system in terms of 
the single particle matrix Green's functions with anoma- 
lous components, 

G(q, ti - t,) = -i{T (^q(ti)^^(i2))>, (3) 

where the four-component spinor 'l'q(t) is 

V'+,. + q(t),VU-q(i))- (4) 

From the equations of motion for the Green's functions 
with Hamiltioninan Ti — fiM, in the Bogliubov approxi- 
mation, where the operators V'-.k ^^^ V'-,k, are replaced 
by \/Nq with Nq the number of condensate particles, and 
with the Hartree-Fock energy included, we obtain 



iq,^J 



(5) 

where G(q.2;) is the Fourier transform of G(q, t) and 

^(q) = q^/^m - fi + g{2no + 2n_ + n+) 
B{q) = {2k + qf/2m - ^ + g{no + n_ + 2n+) 
Z)(q) = 5(-q). (6) 

The chemical potential in leading order is fiQ = 

d{'H)/dNo = grif) + 2gn_ + gn+, where n^ = N^^/V, 
and 
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are the number of particles in the (— ) and (+) states that 
are not in the condensate. In lowest order, 

yl(q) = q^/2m + gn^, B{q) = {2k + qf /2m. (8) 

Low -momentum excitations: The excitation spectra is 
given by the poles of G(q.z) with /i = ^iq; the poles sat- 
isfy 
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Since detG~^(q, 2) = det G~^(— q, — z), the roots 
come in pairs: two positive and two negative, correspond- 
ing to two excitations, for each q. 

One of the two excitations is gapless in the limit q — )• 0, 
and the other is gapless in the limit q — )• —2k. Although 
the roots of det G~^(q, z) = cannot be found analyti- 
cally for general q, we can construct the low energy dis- 
persion relations for momenta around the gapless points. 
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With strong spin-orbit cou- 
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leading order in the low momentum limit | q| ^ k is 
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The dispersion relation for Qy = is linear at low mo- 
menta, as in the usual Bogoliubov spectrum lUSll . Since 
qy/2m is generally larger than 5|An| in typical experi- 
mental setups, the dispersion is essentially quadratic for 
Qx = Qz = 1 26]. On the other hand, the gapless spectrum 
in the limit of small q' = q + 2k is 
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This excitation is quadratic and free-particle like. The 
spectra of the two excitations agree with the result of [19]. 
Condensate depletion: The densities n± are given in 
terms of the Green's functions by 
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where the contour in the z integration surrounds the nega- 
tive poles in the positive sense. To a first approximation we 
neglect n^ inside G(q, z) . We assume that ^J{2mgYnQ 
is small, and in the end we see that An/no is small, justi- 
fying ignoring n± in the integration. 

Before evaluating the integrals in (fT2^ . we show that they 
converge in the ultraviolet. Explicitly, 
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The poles in z at large q behave as -cf /2m, and since at 
large q, A ^ B ^ D ^ q^ /2m naive power counting 
would indicate that the integrand after the z integration be- 
haves as g^^, which combined with three g-integrals yields 
a quadratic ultraviolet divergence. In fact, cancellations in 
the integrand lead to convergence. In evaluating the z inte- 
gral in terms of the poles of the integrand, we employ the 
relation derived from det G^^(q, z) = at the poles. 
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to replace the left side of (fT4ll in the numerator of ([13) with 
the term on the right side. Then 

dzd^'q {gn^f {z-Bf{z + D) 
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In this form the integrand behaves explicitly as q~^ after 
the z-integral, and is thus ultraviolet convergent; the inte- 
gral for n^ is similarly convergent. 

The depletions n^ can be evaluated numerically as a 
function of k/ y/2mgnQ. The number of excited parti- 
cles riex = n_ + n+ is plotted in Fig. [T] Generally, 
n_ ^ n+, and the contribution of n+ to the number of 
excited particles is negligible. As the figure shows, the 
condensate depletion increases with kJ \J2mgnQ, and is of 
order nQy^{2mg)^nQ <^ n, thus justifying our use of the 
Bogoliubov approximation. 
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FIG. 1: The number of excited particles, in units of {2mgnQ)^''^ 
as a function of the spin-orbit coupling strength k in units of 
y^'Imgno- Generally n_ ^ n+ and Uex ~ n^ . 



Ground state energy: The ground state energy density, 
E, in terms of the Green's function is ll20ll 
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where the term in braces is z plus the single-particle 
Hamiltonian in the (— ,+) basis. The integral equals 
gnQ(2mgnQ)'^^^ times a dimensionless function X of 
fi/guQ and K/y/2rngn^. Since the chemical potential 
in mean-field is // = gn^ and n^x is C'((2mgno)'^/^), 
the energy density is /in/2 + Xgn{2mgnY^^. Then, 
writing the chemical potential similarly as /i = gn{\ + 
Y\/{2mgYn), and using /i = dE/dn, one finds Y = 



-lOX; thus the ground state energy is 
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In calculating X, we take fj, = gug; deviations of fj, from 
gno result in higher order corrections. For k — t- 0, one 
finds X = — 1/I5v27r^, which leads to the ground state 
energy derived by Lee and Yang 112 li 12211 . For general k, 
we calculate X numerically. 



Figure |2] shows the shift in the ground state energy, 
AE = E — gn'^/2, in units of {^J {2'mgYn)g'n? /2, as 
a function of K,l^/2mgn^. The energy decreases with in- 
creasing K, and IS.E changes from positive to negative at 
K ~ 0.6-y/2mgno, an effect too small to see in the figure. 
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FIG. 2: The shift in the ground state energy density, Ai5, in units 
of (y^2mg)3n)(;n^/2, as a function of the spin-orbit coupling 
strength k in units of ^2mgnQ. 



EEC at finite temperature: In the absence of interparti- 
cle interactions, bosons with an isotropic Rashba spin-orbit 
coupling do not condense at non-zero temperature because 
the density of states becomes two-dimensional {mK,/2'K) 
at low energy lUql . However, in the presence of interac- 
tions it is possible for the system to Bose condense at finite 
temperature, as we now discuss. 

The number of excited particles at temperature T is 
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where the v sum is over bosonic Matsubara frequencies. 
The system forms a BEG at a given temperature when n^x 
converges in the infrared, and the total particle density ex- 
ceeds Hex- The infrared structure is captured by the z^ = 
component of the Matsubara sum. Since there are two gap- 
less excitations ei(q) and e2(q'): we need to add infrared 
contributions from two limits q — )• and q' — t- 0. In the 
limit of small q and q' , one finds from inverting Eq. (|5}, 
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where C is a constant as q — )• 0. The integral converges in 
the infrared, and thus a BEG can form at finite temperature. 
This result is consistent with Jian and Zhai's effec- 
tive field theory approach to calculate phase fluctuations 



111711 . applied in three dimensions, through the direct re- 
lation between the condensate depletion and the phase 
fluctuations ll23n : 

no ~ n e-«'^W-^('-'))'>/2, \r - r'\ -^ oo. (21) 

Normal state: So far, we have assumed the existence of 
condensate, and proved that the condensate is not destroyed 
by thermal fluctuations. We should also ask whether a non- 
condensed state is favorable at finite temperature. Here 
we obtain the free energy of the normal state within the 
Hartree-Fock approximation and compare the free energies 
with and without a condensate. 

The reduced Hamiltonian within the Hartree-Fock ap- 
proximation with no condensate is 

^HF = -Vg {nt +n\ + n^n+) 

+ E (^^,P V'l.p) (,/^^ 

where A and B are as in (|6]l with Uq = 0. In fact, 
ri- = rij^ = n/2, where n is the total number density 
of particles; namely there is no spontaneous imbalance of 
population in each pseudospin species, as one can prove 
by introducing independent chemical potentials for each 
species, and seeing, as in Appendix B of ll24ll . that the sec- 
ond derivative of the Ginzburg -Landau free energy with re- 
spect to the population imbalance is positive. 
With n± = n/2, the Helmholtz free energy density is 
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where^±(p) = {{p^±Kf +pl} /{2m)- ^ + ^gn/2;i\\e 
chemical potential is determined by the number equation 
n = {l/V) Ep {/(e-(p)) + /(e+(p))}, where f{x) = 
l/{e^^ — 1). Unlike for free bosons in three dimensions, 
it is possible at any temperature to find a value of /i which 
satisfies the number equation, thus the state without con- 
densate is not kinematically forbidden at any non-zero tem- 
perature. 

As T —7- in the absence of a condensate, /i — t- 
3gn/2, and F — t- 3gn^ /4:. This energy is larger than 
the ground state energy with condensate, (gn^ /2){1 + 
0{^y{2mg)^n)). Therefore, at sufficiently low temper- 
ature, a condensate is energetically preferred. At low tem- 
perature F(;U,no) < F(/U,0), so no > 0. The con- 
densate density decreases with temperature, and the tran- 
sition to the normal state, if second order, occurs when 
dF{fi,nQ)/dnQ = at no = 0. Determination of the 
order of the transition, the transition temperature, and pos- 
sible critical exponents at the transition is in progress [12 711 . 
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